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ABSTRACT 


Effects  of  a  model  for  random  bottom  structure  on  acoustic  intensity  in 
isospeed  shallow  water  are  studied.  The  randomness  is  due  to  stochastic 
variations  in  the  bottom  density  and  sound  speed  in  the  horizontal  direc¬ 
tion  beneath  a  plane  water-bottom  interface.  Ray  geometry,  spreading 
loss,  and  bottom  loss  and  phase  shift  are  examined  in  order  to  derive 
formulas  for  mean  intensity  and  the  variance  of  intensity.  The  expres¬ 
sions  obtained  are  sufficiently  general  to  permit  their  use  with  different 
bottom-acoustic  models  of  sound  reflection.  In  this  paper,  for  illustra¬ 
tive  and  comparative  purposes  two  such  models ,  one  developed  by  MacKenzie 
and  the  other  by  Rayleigh,  are  considered.  The  distinctive  acoustic 
consequences  of  bottoms  of  different  density  mean,  variance,  and  horizontal 
correlation  are  discussed,  as  are  comparisons  of  results  for  the  two  bottom- 
reflection  models.  Intensity  moments  are  obtained  also  for  differing  source- 
receiver  range  and  water  depth. 


INTRODUCTION 


The  acoustical  effects  of  random  interface  topography  between  the  water 
and  bottom  has  been  and  is  the  object  of  much  research.  In  contrast,  the 
influence  of  random  variations  in  bottom  structure  has  received  relatively 
less  attention.  One  recent  perspective  on  the  influence  of  these  two  general 
types  of  bottom  effects  on  acoustic  propagation  is  found  in  Ref.  1.  Both 
types  of  bottom  variations  are  particularly  important  in  shallow  water.2 
Indeed,  ocean  regions  may  be  regarded,  for  acoustic  propagation  studies,  as 
shallow  whenever  sound  transmission  is  significantly  influenced  by  the  bottom 
(see  for  instance  Ref.  3).  However,  there  is  such  an  overabundance  of  bottom 
input  parameters  for  intensity  calculations  and  predictions  that  one  model,  or 
even  several,  seem  inadequate  to  account  for  them  all.  For  example,  Rogers4 
lists  twenty-four  separate  inputs  to  a  propagation  loss  model  for  shallow- 
water  acoustics.  Furthermore,  some  of  these  parameters  are  so  poorly  known 
that  it  would  not  be  feasible  to  model  their  acoustical  effects.  Therefore, 
it  has  been  suggested  (for  example,  see  Ref.  4)  that  a  statistical  approach 
could  possibly  give  reasonably  accurate  predictions  for  propagation  loss.  In 
this  paper,  we  use  such  an  approach  to  study  how  horizontal  random  variations 
in  bottom  density  and  sound  speed  affect  sound  transmission. 

The  specific  problem  we  address  assumes  an  ocean  of  shallow  depth,  so 
that  the  water  can  be  assumed  to  have  constant  sound  speed  and  density.  A 
horizontal  water-bottom  interface  is  also  assumed,  and  the  bottom  is  taken  to 
have  small  random  variations  in  sound  speed  and  density  in  the  horizontal 
direction.  Associated  with  these  random  quantities  are  correlations  which 
depend  upon  the  horizontal  distance  separating  two  locations.  The  ray  theory 
of  propagation  is  used  in  the  water,  and  more  than  one  reflection  theory  at 
the  water-bottom  interface  is  studied.  Since  shallow-water  propagation  is 


2 


dominated  by  repeated  bottom  interactions,  we  are  necessarily  considering 
ocean  areas  such  as  the  continental  terrace  region  where  depths  are  less  than 
about  200  Our  primary  objective  is  to  determine  statistics  of  intensity 

at  a  fixed  receiver,  in  terms  of  statistics  of  the  bottom  structure. 

In  Sec.  I  we  derive  expressions  for  per-ray  travel  time,  geometrical 
spreading  loss,  bottom  loss,  and  bottom  phase  shift.  Rather  general  expres¬ 
sions  for  the  mean  intensity  and  its  variance  are  developed  in  Sec.  II,  while 
Sec.  Ill  contains  a  specialization  of  the  results  when  MacKenzie  and  Rayleigh 
bottom-loss  models  are  used.  The  method  by  which  bottom  acoustic  attenuation 
is  incorporated,  and  the  assumed  type  of  horizontal  bottom  correlation  func¬ 
tion,  are  also  discussed.  Section  IV  presents  and  discusses  numerical  results 
obtained  from  the  expressions  we  derived  for  the  case  of  MacKenzie  theory.  In 
Sec.  V,  we  compare  those  results  obtained  previously  using  MacKenzie  theory 
with  those  employing  Rayleigh  theory.  Finally,  Sec.  VI  is  a  summary  of  the 
paper . 

I .  FORMULATION 

We  consider  an  isospeed  ocean  channel  of  constant  depth,  H,  as  suggested 
in  Fig.  1(a).  A  point  sound  source  S  is  located  at  a  depth  hs  and  a  point 
receiver  R  at  depth  hR.  The  horizontal  separation  between  S  and  R  is  denoted 
by  R. 

As  stated  in  the  introduction,  ray  theory  is  used  here,  and  specular 
reflection  at  both  surface  and  bottom  is  assumed.  Rays  are  distinguished  from 
each  other  by  describing  each  by  the  ordered  pair  (n,j).  Thus,  a  ray  is 
denoted  by  rnj ,  where  n  >  0  indicates  the  total  number  of  bottom  reflections 
between  S  and  R.  The  symbol  j  specifies  the  type  of  ray,  of  which  for  n  >  1 
there  are  four,  as  illustrated  in  Fig.  1(b).  For  a  ray  that  leaves  S  in  a 
downward  (or  upward)  direction  and  arrives  at  R  from  below,  j  is  taken  to  be 
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one  (or  two).  Similarly,  j  =  3  and  j  ~  4  correspond  to  the  other  ray-direction 
pairs  at  S  and  R,  as  shown  in  the  figure.  For  the  special  case  n  =  0,  when 
hs  <  hR,  j  takes  on  the  values  3  and  4  but  not  1  and  2;  when  hs  >  hR/  j  can  be 
2  and  4  but  not  1  and  3  [see  Fig.  1(b)]  .  Each  bottom  reflection  of  the  ray  rnj 
can  be  labeled  by  the  unique  triple  (i,j,n),  where  i  numbers  each  bounce 
sequentially  from  one  to  n.  These  three  indices  will  be  essential  in  describ¬ 
ing  two  bottom  parameters  in  our  model.  In  particular,  we  let  and  c^n) 

ij  ij 

denote  the  bottom  density  and  sound  speed  at  bottom  bounce  i  of  ray  rnj .  The 
water  density  p1  and  the  sound  speed  c-j  are  taken  constant.  Returning  to  Fig. 
1(a),  the  angle  that  ray  rnj  makes  with  the  horizontal  as  it  leaves  the  source 
is  denoted  by  0Rj  .  Of  course,  this  is  also  the  magnitude  of  the  inclination 

angle  everywhere  along  r  .  .  We  define  the  inclination  angle  -  JL  <  0  <  JL  to 

nD  2  n3  2 

be  measured  from  the  horizontal,  and  to  be  positive  counterclockwise,  as 
indicated  in  the  figure. 

The  sound  source  5  is  taken  to  emit  a  unity-amplitude  omnidirectional  cw 
signal  sin  cot,  where  to  is  circular  frequency  and  t  is  time.  Then,  ray  rnj 
arrives  at  the  receiver  after  a  time  of  travel  Tnj  in  the  form  Anj  sin(  tot- 0nj  )  . 
The  relative  per-ray  amplitude  Anj  differs  from  unity  largely  because  of 
spreading  loss  and  bottom  loss,  with  other  losses  being  neglected  here.  The 
per-ray  travel  time  and  phase  shifts  at  surface  and  bottom  reflections  cause 
the  phase  <£nj  at  the  receiver  to  differ  from  that  at  the  source.  We  may  write 


nD 


and 


(1  ) 


d>  .  =  iot  . 

n3  nj 


[n  +  1  (A+v)tt]  -Is 


x=1 


(n) 

ij 


(2) 
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In  Eq.  (1),  Anj  is  spreading  loss,  given  by 
A  2 

A  .  =  P  .  '  (3 

nI  nj  ^ 

where  Pnj  is  the  length  of  rnj  .  The  quantity  Bnj  represents  the  decrease  in 
amplitude  due  to  all  the  bottom  reflections  of  ray  rn j ,  and  may  be  written 


n  t  \ 

B  .  =  n  B  n  , 

n3  i-i  ij 

where  is  the  coefficient  of  reflection  of  rn-j 

ij  J 

In  Eq .  (  2 )  , 


(4) 


at  its  ith  bottom  bounce. 


Tnj  pnj/c1  (5) 

is  travel  time,  while  represents  the  phase  shift  at  the  ith  bottom  reflec- 

ij 

tion  of  rnj ♦  The  symbol  X  and  v  are  parameters  defined  by 

X  =  ( 2 . 5- j ) / | 2 . 5“ j  |  (6) 

and 

v  =  (-1)3  , 

which  account  for  surface  reflections  in  the  family  of  rays  with  n 
bounces.  From  geometric  considerations,  a  formula  for  path  length 
be  shown  to  be 

pnj  =  [r2  +  ( 2nH  +  XhR  +  Vhs)2]1/2  . 

Thus,  Anj  is  given  through  Eqs .  (1),  (3),  (4),  and  (8),  and  <j>nj  by 

(5)-(8).  Note  that,  if  the  source  and  receiver  are  located  on  the 
(i.e.,  hR  =  hs  =  0),  then  Eq-  (8)  reduces  to 

Pn1  =  [r2  +  [2nH]2]V2  f  (9) 


bottom 
of  rnj  can 

(8) 

Ecs .  ( 2 )  and 

surface 


so  that,  in  this  particular  case,  the  parameters  X  and  v  are  not  required. 
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Since  the  amplitudes  of  successive  ray  arrivals  converge  to  zero  with 
reasonable  rapidity,  because  of  increasing  bottom  and  spreading  losses,  we  need 
consider  in  total-field  calculations  only  rays  that  experience  some  maximum 
number  N  or  fewer  bottom  reflections.  A  typical  condition  to  determine  N, 
which  we  have  used  in  our  numerical  examples,  is  to  neglect  rays  which  have 
amplitude  less  than  one  percent  of  the  ray  with  largest  amplitude.  Consequent¬ 
ly,  the  total  acoustic  field  at  R,  of  amplitude  A  and  phase  <P,  may  be  written 
in  the  form 

N  4 

A  sin  (u)t-$)  =  \  'i  Anj  sin(u)t-9nj)  .  (10) 

n=0  j  =  1 

It  follows  from  Eq.  (10)  that  A  can  be  written  in  terras  of  the  per-ray  ampli¬ 
tudes  and  phases  as 

N  4  N  4 

A2  =  (  I  l  Anj  sin  <f)nj  )2  +  (  l  l  Anj  cos  <J>nj  )2  .  (11) 

n=0  j=1  n=0  j=1 

We  shall  employ  Eq.  (11)  in  Sec.  II. 

We  turn  now  to  expressions  for  and  S^n^,  where  bottom  properties  at 

ij  ij 

the  location  of  each  ray  reflection  are  taken  to  vary  in  a  small  random  fashion 

about  fixed  "average"  values.  In. our  model  we  consider  variations  both  in  the 

density  p(ny  of  the  bottom  at  the  ith  bounce  of  r  and  in  the  sound  speed  c^n^ 
ij  nj  ij 

of  the  bottom  at  that  point.  Thus,  we  write 

Pij)/Pl  =  ^P2/P1  ^  +  Eij^  '  (12) 

where  p  is  the  constant  horizontally-averaged  density  of  the  bottom,  and 

2  ij 

is  a  small  random  ouantity  [|e(n)|  <<  l).  In  a  similar  manner,  we  let 

ij 


V='"‘  -  (v<=j)(,+5tP 


(n) 


ij 


ID 


(13) 
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where  c  is  the  constant  mean  sound  speed  in  the  bottom  and  | |  <<  1  .  in 
2  ij 

the  next  section,  we  will  discuss  the  first  and  second  moments  of  the  random 
quantities  and  6(n),  including  their  correlation  at  pairs  of  points  along 

ij  ij 

the  water-bottom  interface. 

We  now  expand  the  bottom  loss  in  a  Maclaurin  series  with  respect  to 

ij 

and  6^,  keeping  only  second  and  lower  degree  terms.  It  follows  that 

ij  ij 


8  ( n ) 
ij 


A  .  +  B  .  +  C  .  6  f n  ^  +  d  .  fe(n)}2 

nj  rvj  13  nj  ij  nj  ^  ij  ' 


+  E  ,(e(n)  S(n)]  +  f  , 

nj  v  ij  ij  ;  nj  ^  J 


.  (n)  \2 


(14) 


where  the  coefficients  Anj  through  Fnj  depend  on  the  reflection  theory  and  mean 
bottom  values  that  are  selected.  In  Sec.  IV  we  show  explicit  results  for  the 
reflection  theories  of  MacKenzie  and  Rayleigh.  Of  course,  other  theories  for 
which  these  equations  are  valid  could  also  be  utilized  if  desired.  The  per-ray 
bottom  phase  shift  may  be  expanded  similarly  to  give: 

ij 


s(n)  =  q  +  h  +  I  <5(n)  +  j  (e(n))2 

ij  nj  nj  ij  nj  ij  nj  ij 


(  (n)  An)  \  r  x(n)  \2 

+  K  .  (  e.  .  0.  .j  +  L  .  I  6.  .  J  . 

nj  v  ij  ij  nj  v  11  J 

Again,  the  coefficients  Gnj  through  Lnj  in  Eq.  (15)  are  prescribed  by  the 
chosen  reflection  theory  and  unperturbed  bottom  density  and  sound  speed. 


(15) 


II.  INTENSITY  MOMENTS 

We  turn  now  to  a  consideration  of  the  mean  and  variance  of  the  relative 
intensity  I  at  R.  In  view  of  our  definition  of  A,  I  =  A2  where  A2  is  given  in 
Eq.  (11).  Since  the  per-ray  quantities  Anj  and  <j>nj  in  Eq.  (11)  depend  on  the 
random  variables  and  6^  in  Eqs .  (12)  and  (13),  I  is  necessarily  a  random 

ij  ij 

variable.  Our  first  objective  is  to  express  I  in  terms  of  these  random  variables. 
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We  begin  by  expressing  <j>nj  from  Eq.  (2)  in  terms  of  quanitities  <J>n^  and 


nj 


Y  ,  where  <£>  does  not  depend  upon  the  small  stochastic  perturbations 
nj  nj  ij 

and  while  y  does: 

ij  nj 


^nj  ~  ^nj  +  ^nj 


(16) 


To  second-degree  terms,  Eqs .  (2),  (5),  (8),  and  (15)  give 


$  .  =  wp  ./c  -  [n  +  1.  (A+vIItt  -  n  G  . 

nj  nD  1  9  v  J  J  uj 


(17a) 


and 


Tnj 


l  [h  .  £(n)  +  I  .  6(n)  +  j  .  0(n)  ]2  +  K  .  (e(n)  6(n)  ) 

i^i  n]  nj  ID  nD  v  id  J  id  id 


+  L  .(6(n))2] 

nD  v  id  J 


(17b) 


Next,  we  write  Anj  as 


^nj  ^nj  ^nj 


(18) 


where  similarly  A  is  independent  of  and  and  T  is  one  when  the 

nj  ij  ij  nj 

small  perturbations  are  zero.  From  Eqs,  (1),  (3),  (4),  and  (14),  we  find 


nD 


An./P  . 
nD  nD 


(19a) 


and 


r  . 

nD 


=  n  [l  +  (b  ./A  .)efn)  +  (c  ./A  .  )6(r 


,(n) 


i=1 


nD  nD  id 


nD  nD  '  ID 


+  (D  /A  +  (e  /A  )(e^  6^) 

nj  nj  ij  n~i  ni  ii  ii  ' 


nD  nD  id  id 


+  (?  ./A  0(«<n>)2]  ■ 

v  nj '  nj  '  v  l  j  J  J 


(19b) 
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Substituting  Eqs .  (16)  and  (18)  into  Eq.  (11),  carrying  out  the  squares, 
and  combining  terms,  we  obtain: 


N  N  4  4  __ 

i  =  i  i  i  r  a  . 

n=0  m=0  j=1  k=1  n3 

With  an  expression  of  the  cosine  that  retains  all  terms  to  second  order  in  the 
density  and  sound-speed  perturbations,  it  follows  from  Eq.  (20)  that 

N  N  4  4  _ 

i  =  l  l  l  l  a  . 

n=0  m=0  j=1  k=1  n3 


A  r  .  r  cos 

mk  nj  mk 


(*  . ,  +Y  .  -T  ,  )  • 

mk  nn  mk  J 


(20) 


x  cos  ($  . -<&  )  +  (y  .-Y  .  )r  .  T  sinf  $  . ,  )}  .  (21 

n  j  mk J  nj  mk J  n j  mk  ^  nj  mk  J 1 

Thus,  the  intensity  is  known  from  Eqs.  (8),  (17),  (19),  and  (21),  correct  to 

second-order  terms  in  and  6^n^. 

ij  ij 

We  are  now  in  a  position  to  determine  the  mean  and  variance  of  the 
intensity, 


y(i)  =  e(i)  (22) 

and 

a2(l)  =  E( I2 )  -  E2(I)  ,  (23) 

where  E  denotes  stochastic  expectation.  In  doing  so,  we  will  take 

E(£(n))  =  E(6(n)J  =  0  ,  (24) 

V  XJ  J  K  XJ  ' 

o2(e(n))  =  a2  ,  (25a) 

ij  e 


and 


o 2{6(n)) 


(25b) 
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Equations  (25)  mean  that  all  and  all  are  assumed  to  have  the  same 

ij  ij 

variance.  Other  second-order  moments  will  be  written  as 


and 


.(n) 

£0(m))  = 

a2 

ij 

e 

■  (n) 

6(m)l  = 

2 
0  o 

*£)  - 

’ij 

Ik  } 

6 

;  ( n ) 

II 

e 

0 

a  C(e(n>,  6(a) 

ij 

£k 

£ 

6  ij  £k 

(26) 

0  and 
e 

a 

6 

are  the  standard 

(26a) 


(26b) 


(26c) 
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duced  in  Eqs .  (25),  and  C(a,b)  represents  the  correlation  coefficient  of  the 
random  quantities  a  and  b. 

In  calculating  moments  involving  the  we  employ  for  convenience  an 

ij 

approximate  relationship  between  bottom  sound-speed  and  density  perturbations. 
For  small  enough  perturbations,  it  may  be  anticipated  that  and  are 

ij  ij 

approximately  proportional,  which  is  confirmed  by  examination  of  data.  In 
particular,  from  the  extensive  data  in  Table  I  of  Ref.  6  and  linear  regression, 
or  using  appropriate  equations  from  that  reference,  we  found  that  to  a  good 
approximation, 


£(n)  =  -  20(n) 
ij  *D 


(27) 


This  result  can  be  used  to  write  Eqs.  (26b)  and  (26c)  as 
2  ^  r  r  (n)  r(m)^  i  2  Jn) 


ID 


p  {  x  l  n  i  -(m)1!  i  z  p  f  in  j 


ID 


.  (m) 

'£k 


(28a) 


and 


a  Ojt 

e  6 


c(« 


(n) 


= 

l  k  ; 


- 1  a2  c(£(n),  e;m)) , 
2  e  i]  ikJ 


(28b) 
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so  that  only  the  single  correlation  function  C(e^n^,  )  enters  our  calcula- 

ij  £k 

tions . 

We  can  now  perform  the  stochastic  averaging  and  write  expressions  for  p( I) 
and  a2(l)  from  Eqs.  (22)  and  (23),  by  use  of  Eqs.  (17b),  (19b),  (21),  (24)- 

(26),  and  (28).  These  formulas,  while  important  to  our  analysis,  are  quite 
lengthy,  and  thus  will  not  be  presented  here.  Instead,  we  next  further  simplify 
the  mean  and  variance  of  the  intensity,  which  after  the  stochastic  averaging 
still  contain  the  functions  cos ( $nj -S^)  and  sin(  •  We  do  this  by  per¬ 

forming  an  average  over  those  phase  contributions  which  are  independent  of  our 
stochastic  bottom  variations.  This  procedure  is  essentially  equivalent  to  the 
process  of  obtaining  incoherent  total-field  intensity.  Specifically,  we  assume 
that  the  $nj ,  and  differences  between  them,  are  random  and  uniformly  distributed 
over  a  2ir-radian  interval.  Performing  such  phase  averages  and  denoting  them  by 
angular  brackets,  we  obtain 
N  4 

<P(I)>  =  l  l  A2.  (29) 

n-0  j  =1  nD 

and 


Ol 

mk 


i=1 


l  c n  , 

ID 


£= 1 


(30) 


It  should  be  recalled  that  in  Eqs.  (29)  and  (30),  n  or  m  equal  to  zero  is  a 
special  case  with  only  two  possible  j  values  (see  discussion  at  start  of  Sec.  I) . 
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Equations  (29)  and  (30)  represent  the  principal  results  of  this  section 
and  will  be  analyzed  in  subsequent  sections.  At  this  point  we  content  our¬ 
selves  with  two  important  observations.  First,  the  mean  intensity  <p(I)> 
depends  on  mean  per-ray  amplitudes  only,  and  does  not  vary  with  any  of  the 
bottom  coefficients  Bnj , . . . , Fnj , Hnj , . . . ,Lnj  from  Eqs .  (14)  and  (15).  Second, 

the  variance  <o2(I)>  depends  only  on  the  coefficients  Bnj  ,  Cn j  ,  Hn j ,  and  Inj  of 
linear  terms  in  and  6^  in  Eqs.  (14)  and  (15).  That  is,  even  though  the 

ij  ij 

variance  is  a  second-order  quantity  in  the  density  and  sound-speed  perturba¬ 
tions,  and  even  though  six  coefficients  of  second-order  terms  (i.e.,  Dnj ,  En j , 
Fn j  r  Jn j '  Kn j '  and  Lnj ^  are  retained  in  calculating  Eqs.  (29)  and  (30),  our 
final  results  are  independent  of  these  coef f icients .  This  is  a  significant 
simplification,  since  these  six  coefficients  represent  certain  higher-order 
information  about  the  acoustic  response  of  the  bottom,  which  is  thus  not 
required  for  our  intensity  results . 

III.  BOTTOM  MODELING 

Having  derived  formulas  for  the  mean  intensity  <u(I)>  and  variance 
<a2(I)>,  stochastically  averaged  over  random  bottom  variations  and  incoherently 
averaged  over  per-ray  phase,  we  wish  to  use  them  with  two  models  of  bottom 
reflection . 

We  consider  first  MacKenzie's  model,7  in  which  bottom  loss  is 


(31a) 


and  bottom  phase  shift  is 


tan  1  {2ah  sin  0  .  [a2  sin  c  .  -  fh2  +  h2 1 ]  1  }  . 

1  nj  1  nj  K  J 1  J 


(31b) 


In  Eqs .  (31) 
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h+  =  h  9  +  [(<*/&) 2  +  g2]1/2}1/2  , 

2.  2.  p 

g  =1  [l  -  (v  /c  )cos  6  -  (a/g)  ]  , 

2  z  1  n3 


(32a) 

(32b) 


cos  6  .  =  R/P  . 
n:  nj 


(32c) 


and 


a 


Vpici 


( 32d) 


The  quantity  V2  is  the  phase  speed  of  acoustic  waves  in  the  bottom  and,  to  an 

excellent  approximation, 7  v  =  c^n)  .  The  term  a/3  is  a  factor  that  accounts 

2  ij 

for  attenuation  in  the  bottom,  and  in  fact,  a  is  often  referred  to  as  the 
attenuation  coefficient.  This  quantity  provides  a  dissipative  mechanism  for 
sound  which  enters  the  bottom  and  is  typically  a  function  of  bottom  density. 
However,  since  the  values  of  a/3  'are  usually  extremely  small^  compared  to  other 
quantities  in  our  formulas  for  <jj(l)>  and  <a2(I)>,  we  may  reasonably  express 
a/3  in  terms  of  the  mean  density  p2,  rather  than  in  terms  of  .  That  is, 

ij 

the  perturbation  factor  in  is  negligible  in  a/3.  We  now  turn  to  a 

ij  ij 

model  for  the  dependence  of  this  quantity  on  p2* 

Hamilton^  gives  an  expression  kp(D/2Tr  for  a  in  terms  of  a  quantity  kp, 
which  characterizes  natural  saturated  surface  sediments,  and  the  circular 
acoustic  frequency  ca.  His  findings  are  presented  in  the  form  of  scatter 
diagrams  and  regression  curves  for  kp  versus  porosity.  Since  porosity  can  be 
expressed  fairly  accurately  in  terms  of  bottom  density, ^  we  can  obtain 
expressions  for  kp  in  terms  of  P2*  We  found  from  Hamilton's  results  that: 

kp  =  1  .1  096  -  0.3058  p2  ,  (33a) 

for  coarse,  medium,  and  fine  sand; 


kp  =  5.9990  -  2.845  p2  , 


(33b) 
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for  very  fine  sand  and  lower  porosity  mixed-size  sand; 

kp  =  -4.4240  +  2.837  P2  ,  (33c) 

for  mixed-size  sand;  and 

k  =  0.3622  -  0.5712  p„  +  0.2626  p2  *  (33d) 

p  2  2 


for  silt-clay  bottoms.  The  quantity  g  is  given  by^ 
our  earlier  approximation,,  we  may  write 


U)/c(n)  / 
ij 


and  consistent  with 


a/0  =  k  /  2tt  =  k  c  /2tt  .  (34) 

P  i]  P  2/ 

Note  that  a/g  from  Eqs .  (33)  and  (34)  is  independent  of  frequency.  Other 
investigators  have  proposed  different  models  for  attenuation, 1 1 ^  and  they 
could  also  be  used  if  desired. 

From  our  bottom  model,  we  must  obtain  the  coefficients  Anj  ,  -  Cnj/2, 

and  Hnj  -  Inj/2  for  use  in  Eq.  (30).  It  follows  from  Eqs.  (31)  and  (32)  and 
our  assumption  on  a/g  [see  Eqs.  (33)  and  (34)]  that  the  variations  and 

ij 

6(n)  in  p(R)  and  c(R)  only  enter  Eqs.  (31)  and  (32)  through  the  quantities 

ij  ij  ij 

c(n)/c  in  £q.  (32b)  and  the  a  in  Eq.  (32d).  This  latter  parameter  is  some- 

ij  1 

times  referred  to  as  the  impedance  ratio,  i.e.  the  ratio  of  acoustic  impedance 
of  the  bottom  to  that  of  the  water.  Thus,  we  write  c  /c(R)  and  p(R)/p  as  in 

1  ij  ij  1 

Eqs.  (12)  and  (13),  and  with  Eq.  (27),  we  can  express  /c  and  a  in  terms 

ij  1 

of  5(n)  and  the  mean  sound-speed  and  impedance  ratios 

ij 

C  =  <=x/c,  (35a) 

and 


a  =  P2C2/P1C1  .  (35b) 

In  addition  we  relate  the  bottom  sound  speed  c2  (in  m/sl  in  our  model  to  p2  (in 
g/cm3)  by  the  relation6 
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c2(p2)  =  2270.9  -  1194.4  P2  +  474.6  . 


(35c) 


We  then  expand  Eqs .  (31a),  (32a),  and  (32b)  in  a  Maclaurin  series  in  6^. 

ij 

With  the  notation  of  Eq.  (14),  we  find 


A  .  =  M  ./M  .  ' 

nj  nj  nj 


and 


C  .  -  2B  .  =  )  -  Q  .  M  ./(m+.)3]/2  , 

nj  ltvj  K  nj  113  '  nj  113  v  n]  J  J 


(36a) 


(36b) 


where 


II 

+1 

£ 

[2T  .  ± 

20  R+. 

sin  0  . 

-2  .  2 
+  a  sm 

6  .] 

nj 

n: 

nD 

nj  J 

O 

i+ 

II 

t-. 2  2 

[c  cos 

e  .)( 

S  . /T  .  ) 

—2  2 
-  6a  sin 

0  . 

nj 

n  j  J  ' 

n:  nj 

n3 

1/2 


(36c) 


—  + 


+  a  R  .  sin  6  .  ffc2cos20  .  )/T  .  ±  6*1  , 
nu  nj  L  ^  nj  J  nj  J 


± 

R  . 
nD 


n: 


It  .  ±  s  ./2] 

L  nj  nj'  J 


1/2 

/ 

2  .  _2  , .  1 1/2 


[(o/e)  +  si/4]  7  , 


( 36d ) 

( 36e ) 

( 36f ) 


and 


S  .  =  1  -  c2cos0  .  -  fa/ftV 
nj  nj  K  J 


(  36g ) 


These  coefficients,  though  lengthy,  present  no  difficulty  for  numerical  evalua¬ 
tions.  In  a  similar  way,  we  find  from  Eqs.  (31b),  (32a),  and  (32b)  that 

I  .  -  2H  .  =  [(-X  .U  ,+V  .W  .  )/fb 2.+V\i)]  ,  (37a) 

nj  n3  nj  nj  nj  n]  ^  nj  n;j  J 


U  .  =  c2sin20  .  -  2T  .  , 
n3  n]  ,n] 


V  =  c^cos^0  [s  /T  )  +  6a2sin^0  , 
nj  nj  nj  nj  nj 


(37b) 

(37c) 


and 


W  .  =  20  R  .sin0  .  , 

nj  n]  nj 


_  2  ^ 

X  .  =  a  R  . sin0  . [ ( c  cos  0  .  )/T  .  +  6  ]  . 
n]  n]  n;jLV  n];  n]  J 


(37d) 


( 37e ) 
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A  second  and  simpler  bottom-reflection  model  is  that  obtained  from  the 
classical  theory  of  Rayleigh.12  Indeed,  Rayleigh  bottom  loss  and  phase  shift 
can  be  found  from  the  corresponding  MacKenzie  results,  Eqs .  (31)  and  (32),  by 
setting  a/3  =  0-  In  this  case,  the  bottom  loss  can  be  written  as 


b  : 


( n) 


"  l(pln,/pj  -  [(c,/c'”T  -  <=os"e  -  cos2e 

1  v  in  1  1  J  L  v  1 '  in  J  nn  J  1  m  J  1 


(n)  >2 


1/2  | 


1-1/2 


x  ((p'f/p,)  +  [(=,/«!"’)"  -  cos^e  ]V^[,  _  cos2 e  r/2r’  - 


(n)  -\2 


1/2 


1-1/2  1-1 


1=) 


&c  <  e  <  tt/2  , 


( 38a) 


and 


B  {n) 

ID 


=  1  , 


0  <  0  <  0 


(38b) 


For  the  bottom  phase  shift,  we  may  write 


s(n] 

ID 


=  2tan  1  { I cos^Q  . 

1 L  1  ni 


n: 


-  [v“)!],/![»2('-'“»nj),"r,/'l.  0  <  8<  9c 


2.  ^1/2 ,-1/2 


(39a) 


and 


S(n)  =  0  ,  0  <  0  <  TT/2  .  (39b) 

ij  c 


In  Eqs.  (38)  and  (39),  the  critical  angle  is  defined  by1 2 


cosQ  =  c  /c|n^  ,  (40) 

c  1  13 


and  0nj  is  given  by  Eq.  (32c). 

Once  again,  by  writing  c  /c^n^  and  p^/P  as  in  Eqs.  (12)  and  (13)  and 

1  ij  ij  1 

using  Eq.  (27),  we  can  expand  Eqs.  (38)-(40)  in  a  Maclaurin  series  in  as 

ij 

done  for  the  MacKenzie  model.  This  procedure  determines  the  quantities 
(Bnj-Cnj/2)/Anj  and  Hnj-Inj/2  which  appear  in  Eq.  (30).  We  do  not  display 


these  results  because  they  can  be  obtained  from  the  case  a/3  =  0  of  Eqs.  (36) 


and  (37 )  . 
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|  n 

For  evaluation  of  the  intensity  variance  <az(l)>  in  Eq.  (30),  we  also  need 
the  correlation  coefficient  in  Eq.  (28).  For  convenience  we  will  use  the 
Gaussian  form^ 

n(  (n)  (m) \  r  r  r  (n)  (m)^  ]2i 

C(eij  '  kl'*  =  eXpt"  Xij  "  ^  '  '  (41) 

where  (x^n^  -  x^m^)  is  the  horizontal  distance  between  the  ith  bottom  reflec- 
kJL 

tion  of  ray  rnj  and  the  kth  bottom  bounce  of  ray  rm£.  The  parameter  L  is 
called  the  correlation  length,  which  is  the  value  of  |x^n^  -  x^m^ |  where  C 

ij  k£ 

takes  the  value  of  1 /e .  Of  course,  correlation  functions  other  than  that  in 
Eq.  (41)  could  be  used,  if  desired  or  if  suggested  by  bottom  data. 

IV.  MACKENZIE- BOTTOM  RESULTS 

We  now  discuss  some  results  obtained  using  our  expressions,  Eqs .  (29)  and 
(30),  for  mean  intensity,  <p(l)>,  and  intensity  variance,  <o^(I)>.  Equations 
(33)  and  (34)  were  used  to  determine  values  for  a/3,  Eqs.  (36)  and  (37)  for 
coefficients  in  Eq.  (30),  and  the  correlation  coefficient  in  Eq.  (41)  was 
employed.  For  convenience,  our  numerical  results  were  determined  with  the 
source  and  receiver  located  on  the  surface.  .  Of  course,  other  configurations 
could  have  been  used.  In  all  our  computations,  we  found  that  all  rays  were 
included  up  to  a  maximum  number  N  of  6  or  7  bottom  bounces,  according  to  the 
criterion  stated  in  Sec.  I.  Rays  with  more  bottom  interactions  were  shown  to 
negligibly  affect  our  results. 

Graphs  of  the  function 

Mr  =  10  log o  <U  (I)  >/<a(l)  >  (42) 

appear  in  Figs.  2  and  3  and  will  be  discussed  shortly.  In  addition,  graphs  of 


intensity  ratios  IR/  where  this  symbol  is  used  to  denote  both 
IR  =  10  log -j  o  <  u  ( I)  >/<U0(  I)  > 


(43a) 
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and 

IR  =  10  log-j  o (  <u  ( I)  >  ±  <cr ( I ) >  ) /< Hq ( I ) >  ,  (43b) 

are  considered  in  Figs.  4-7.  The  moment  ratio  MR  is  the  ratio  in  dB  of  mean 
intensity  to  the  standard  deviation  of  intensity.  It  can  be  thought  of  as  a 
signal-to-noise  ratio,  in  the  sense  that  a  small  (i.e.,  a  large  negative)  value 
of  Mr  means  that  the  fluctuation  in  intensity  is  large  relative  to  the  mean;  a 
large  value  corresponds  to  the  standard  deviation  being  small  compared  to  the 
mean.  Equation  (43a)  is  a  normalized  dB-measure  of  intensity,  in  which  the 
term  <y0(l)>  is  taken  to  be  the  mean  intensity  for  a  bottom  density  with  the 
particular  value  p2  =  1.80g/cm3.  This  value  corresponds  to  a  typical  silty- 
sand  bottom.  In  Eq.  (43b),  < y q ( I ) >  has  the  same  interpretation  as  in  Eq. 

(43a);  Eq.  (43b)  itself  is  a  normalized  dB-measure  of  intensity  spread.  That 
is,  the  two  expressions  in  Eq.  (43b)  describe  the  amount  of  intensity  varia¬ 
tion  within  one  standard  deviation  of  the  mean,  where  the  mean  is  given  by  Eq- 
(43a).  A  small  difference  in  the  values  of  the  two  quantities  in  Eq.  (43b) 
indicates  that  the  fluctuations  in  intensity  are  small,  while  a  large  differ¬ 
ence  suggests  a  significant  fluctuation. 

Figure  2  is  a  plot  of  the  moment  ratio,  Eq.  (42),  versus  mean  bottom 

density  p2*  The  values  used  for  p2  are  typical  of  a  bottom  on  the  continental 

shelf.  Density  and  sound-speed  values  in  the  water  were  taken  to  be  p-j  = 

1.025  g/cm3  and  c1  =  1523  m/s.  In  this  figure,  as  in  all  our  graphs  except 

Fig.  7,  we  have  chosen  oz  =  0.1.  This  relation  was  made  from  an  examination  of 

Table  I  of  Ref.  6.  However,  it  should  be  pointed  out  that,  due  to  the  manner 

in  which  a2  appears  in  Eq.  (30),  graphs  of  Mo  would  merely  be  shifted  verti- 
e 

cally  if  ae  were  assigned  a  different  value.  Thus,  the  shapes  and  relative 
values  of  the  curves  would  not  be  altered.  In  Fig.  2,  we  have  taken  L/H  =  1, 
so  that  the  bottom  correlation  length  has  been  selected  to  be  equal  to  the 
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water  depth.  Other  numerical  calculations,  some  of  which  are  discussed  in 
connection  with  a  later  figure,  have  explored  the  relatively  small  influence  of 
different  values  for  L/H. 

Figure  2,  then,  illustrates  the  behavior  of  mean  intensity  as  a  function 
of  mean  bottom  density,  for  four  values  of  range-to-depth  ratio  R/H.  We 
observe  that  MR  does  not  appear  to  exhibit  a  well-defined  behavior  as  R/H 
changes.  However,  values  of  the  moment  ratio  can  be  substantially  different 
for  different  aspect  ratios  R/H.  For  example,  when  p2  =  1.6/cm3,  the  MR  value 
for  R/H  =  10  exceeds  that  for  R/H  =  5  by  about  8  dB.  All  MR  versus  p2  curves 
exhibit  a  general  upward  trend,  indicating  that  the  mean  intensity  tends  to 
increase  as  mean  bottom  density  increases.  As  an  illustration,  MR  increases  by 
nearly  20  dB  as  P2  increases  from  1.55  to  2.10  g/cm3.  As  described  earlier, 
this  behavior  suggests  a  much  greater  fluctuation  in  received  intensity  for  a 
low-density  bottom,  such  as  one  composed  of  silty  clay.  Rogers4  obtained  the 
same  anomalous  type  of  conclusion  using  an  algebraic  model  for  propagation 
loss.  Specifically,  his  loss  model  is  based  on  a  combination  of  information 
from  Weston,14  Urick,15  and  McPherson  and  Daintith.1^  However,  Ref.  4  utilizes 
a  combination  of  ray  and  mode  theories  of  propagation  over  a  deterministic 
bottom,  while  we  employ  ray  propagation  over  a  horizontally  random  bottom. 
Returning  to  Fig.  2,  one  way  to  explain  the  increase  in  MR  with  p2  is  to  con¬ 
sider  properties  of  Mackenzie  reflection  theory.7  In  particular,  the  change  in 
bottom  loss  and  bottom  phase  shift  with  increasing  grazing  angle  is 

ij  ,  ij 

much  greater  for  slow  (low  density)  bottoms  than  for  fast  (high  density) 
bottoms.  This  is  particularly  true  for  small  ray  angles,  which  are  dominant  in 
isospeed  channels  with  large  aspect  ratios.  Thus,  these  greater  changes  would 
account  for  larger  values  of  <a(I)>  relative  to  values  of  <p(l)>  when  p2  is 


small . 
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In  Fig.  3  we  have  kept  the  aspect  ratio  constant  (R/H  =  20)  and  varied 
L/H,  in  graphs  of  moment  ratio  versus  mean  bottom  density.  It  is  readily 
observed  that  as  p2  increases,  MR  exhibits  an  overall  increasing  trend.  When 
the  normalized  correlation  length  is  taken  as  L/H  =  1,  for  example,  MR  is 
observed  to  increase  by  about  14  dB  as  p2  increases  from  1.42  to  2.10  g/cm3. 

We  also  find  little  dependence  of  MR  on  L/H,  although  the  moment  ratio 
tends  to  be  larger  for  smaller  L/H  values  when  p2  is  large.  However,  MR  does 
not  behave  monotonically  with  L/H.  As  L  becomes  progressively  larger,  the 
correlation  coefficient,  Eq.  (41),  increases.  Thus,  for  a  small  value  of  L, 
the  correlation  between  random  bottom  properties  at  two  points  is  small,  and 
the  fluctuations  in  the  densities  at  these  points  tend  to  cancel.  Conse- 

ij 

quently,  <a(l)>  decreases,  causing  MR  to  increase  with  L  when  L  is  small,  as 
suggested  by  the  L/H  =  0  and  1  curves  in  Fig.  3.  On  the  other  hand,  for  larger 
L,  there  is  more  of  a  density  relationship  at  two  bottom  points.  The  reduc¬ 
tion  in  fluctuations  will  not  occur  to  as  great  an  extent  as  for  small  L,  so 
that  the  value  of  <a(I)>  tends  to  be  larger  for  larger  L.  This  results  in  a 
decrease  in  MR  with  increasing  L,  as  observed  in  the  L/H  =  10  and  00  curves  of 
Fig.  3.  Indeed,  in  our  model,  L/H  =•  00  corresponds  to  a  horizontally-unif orm 
bottom  of  density  which  is  unknown  but  is  close  to  p2.  Finally,  we  note  that 
calculations  with  R/H  values  other  than  20  were  found  to  have  little  gross 
effect  on  the  behavior  of  MR  with  p2  for  the  selected  L/H  values. 

Figures  4  through  7  are  graphs  of  intensity  ratios,  as  given  in  Eqs .  (43). 
As  we  discussed  just  following  that  equation,  these  ratios  illustrate  the 
amount  of  spread  in  intensity  predicted  by  our  model  for  different  values  of 
p2,  R/H,  and  L/H.  The  quantity  < y 0 ( I )  >  is  the  mean  intensity  for  mean  bottom 
density  p2  =  1.80g/cm3,  and  is  used  as  the  reference  value  in  the  dB  plots. 
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In  Figs.  4  and  5,  the  normalized  bottom  correlation  length  is  unity,  but 
the  aspect  ratio  R/H  is  5  in  the  former  figure  and  20  in  the  latter.  The 
dashed  curves  describe  mean  intensity,  and  the  vertical  separation  between 
solid  curves  represents  the  amount  of  intensity  fluctuation  within  one  standard 
deviation  about  the  mean.  In  both  figures,  the  dashed  intensity  curves  show  a 
general  upward  trend  as  p2  increases,  indicating  that  mean  intensity  increases 
as  the  bottom  becomes  increasingly  fast.  Both  figures  show  that  intensity 
variation  is  much  greater  for  lower  density  bottoms  than  for  higher  density 
ones.  This  is  consistent  with  our  previous  discussion  of  Fig.  2.  In  Fig.  4, 
for  example,  the  spread  between  the  <p(I)>  ±  <o(l)>  curves  is  nearly  10  dB  when 
p2  =  1.70.  The  spread  in  values  of  intensity  (solid  curves)  tends  to  be  great¬ 
er  in  Fig.  5  than  in  Fig.  4,  so  that  intensity  variation  appears  to  increase 
with  R/H. 

Figure  6  displays  the  effect  on  intensity  ratios  of  the  normalized  corre¬ 
lation  length  L/H.  Increasing  the  length  is  seen  to  widen  the  dB  interval  in 
which  intensity  is  within  one  standard  deviation  of  its  mean.  Again,  we 
attribute  this  to  the  cancelling  of  fluctuations  in  bottom  density,  as  dis¬ 
cussed  in  conjunction  with  Fig.  3.  For  example,  when  p2  =  1.90,  the  spread  is 
3.5  dB  as  L/H  approaches  zero,  but  is  about  6.0  dB  as  L/H  becomes  very  large. 
Changing  the  value  of  L  appears  to  have  a  slightly  greater  effect  for  trans¬ 
mission  over  higher  density  bottoms.  In  slow  bottoms  the  effects  of  large 
changes  in  bottom  loss  and  phase  shift  with  grazing  angle,  discussed 

ij  ij 

previously,  causes  the  large  interval  between  the  <y>  ±  < a>  curves. 

The  effect  on  intensity  variation  of  changing  a£,  the  standard  deviation 
of  bottom-density  variations  in  Eq.  (25a),  is  illustrated  in  Fig.  7.  Here,  L/H 
and  R/H  are  held  constant,  and  a£  is  assigned  four  different  values.  Increasing 
ae  is  observed  to  cause  widening  of  the  spread  of  the  intensity  distribution. 
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This  is  reasonable  since  <o>  increases  with  ae,  causing  a  greater  variance  of 
the  signal  intensity.  As  an  example,  the  one-standard-deviation  spread 
increases  from  0  dB  when  a£  =  0  to  about  13  dB  when  o£  =  0.20,  for  p2  = 
1.80g/cm3.  The  zero  intensity  value  corresponding  to  o£  =  0  results,  of  course, 
from  the  fact  that  there  is  no  intensity  fluctuation  when  the  bottom  has  con¬ 
stant  density.  Finally,  we  note  that  the  effect  of  increased  oe  is  greatest 
when  dealing  with  low  density  bottoms. 

V.  RAYLEIGH- BOTTOM  RESULTS 

Up  to  this  point,  we  have  used  MacKenzie  reflection  theory  in  examining  the 
effects  of  bottom-structure  variations  on  received  acoustic  intensity.  In  Figs. 
8  and  9,  however,  we  compare  intensity  results  using  both  MacKenzie  and  Rayleigh 
bottom  models.  The  former  employs  our  expansions  Eqs.  (36)  and  (37)  of  Eqs . 

(31)  and  (32),  while  the  latter  uses  analogous  expansions  of  Eqs.  (38)  and  (39). 
In  Fig.  8  we  plot  the  moment  ratio  MR,  Eq.  (42),  versus  p2  for  R/H  =  5  (solid 
curves)  and  20  (dashed  curves),  while  fixing  L/H  at  the  value  one.  The  heavy 
curves,  which  appeared  previously  in  Fig.  2,  are  the  results  for  the  MacKenzie 
bottom  model.  The  light  curves  correspond  to  the  Rayleigh- theory  bottom.  In 
Fig.  9  we  plot  MR  versus  p2  for  L/H  =  0  (solid)  and  L/H  +  «  (dashed),  when  R/H 
is  fixed. 

We  observe  that  the  light-dashed  curves  in  Figs.  8  and  9,  and  the  light- 
solid  curve  in  Fig.  9,  approach  vertical  asymptotes  at  several  values  of  p2. 

The  light-solid  curve  in  Fig.  8  has  a  single  vertical  asymptote.  In  contrast, 
moment-ratio  curves  for  a  MacKenzie  bottom  do  not  exhibit  this  behavior.  A 
brief  explanation  of  these  asymptotes  is  appropriate  here.  It  is  known  that  for 
c2  >  c-j  and  p2  >  p-j,  the  Rayleigh  bottom  reflection  coefficient  increases  to  a 
value  of  unity  as  the  angle  of  incidence  varies  from  normal  to  the  critical 
angle  0C.  At  angles  more  grazing  than  0G,  the  relative  amplitude  of  the 
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reflected  wave  remains  one,  but  there  is  a  phase  change  between  the  reflected 
and  incident  waves,  given  by  Eqs.  (40). 12  It  follows  from  Eqs *  (32c)  and  (40) 
that 

[(R/Pnj)2  -  (ci/c2)2]  =  0  (44a) 

when  the  critical  angle  0C  is  attained.  Equation  (44a)  is  equivalent  to 

c2  =  cl0  +  ( 2nH/R.}  2]  1/2  #  (44b) 

where  we  have  used  the  path-length  formula,  Eq.  (9).  In  our  expansions  of 

and  S(n>,  Eqs.  (38)  and  (39),  the  left  side  of  Eq.  (44a)  occurs  in  the  * 

ij  ij 

denominators  of  the  terms  Cnj  and  IRj .  Thus,  for  given  number  n  of  bottom 
reflections  and  given  range  R,  a  value  of  C2  (and  hence  of  p2,  since  they  are 
related  by  Eq.  (35c))  can  occur  for  which  Cnj  and  IRj  become  infinite.  This 
feature  is  also  apparent  in  the  graphs  of  B^  and  where  it  is  manifested 

ij  ij 

as  the  infinite  slopes  at  the  critical  angle.12  Therefore,  asymptotes  of  the 
M2  versus  p2  curves  in  Figs.  8  and  9  result  from  our  expansions  in  and 

ij 

6^),  which  are  not  valid  at  the  critical  angle.  Because  C  and  I  become 
ij  nj  nj 

infinite,  so  does  <a>,  so  that  MR  expressed  in  dB  approaches  negative  infinity. 
If  the  exact  (and  more  complicated)  dependence  on  and  6^n^,  rather  than 

ij  ij 

our  linear  approximations,  were  in  some  asymptotic  manner  retained  in 

ij 

and  S<n) ,  then  the  moment-ratio  results  would  not  be  expected  to  possess  dis- 

ij 

continuities  in  values  and  slopes.  Instead,  it  would  be  anticipated  that  MR  in 
dB  would  still  exhibit  sharp  decreases  at  certain  p2  values,  but  these 
decreases  would  not  become  infinite.  We  remark  also  that  as  R/H  increases, 
more  asymptotes  (corresponding  to  different  n  values)  appear  in  the  MR  versus 
P2  curve,  as  implied  by  Eq.  (44b). 

Many  of  the  same  features  pointed  out  in  the  discussion  of  Fig.  2  can  be 
seen  in  Fig.  8.  For  instance,  we  observe  that  MR  exhibits  a  rather  complicated 
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behavior  as  R/H  changes,  for  either  MacKenzie  or  Rayleigh  theory.  However,  the 
values  of  the  moment  ratio  can  be  quite  different  for  different  R/H  values. 

For  example,  when  p2  =  1.85,  the  Rayleigh  MR  value  for  R/H  =  5  exceeds  that  for 
R/H  =  20  by  about  8  dB.  The  difference  at  the  same  p2  value  using  the 
MacKenzie  curves  is  about  5  dB.  Apart  from  the  previously-noted  discontinui¬ 
ties  in  the  Rayleigh  curves,  there  are  some  general  qualitative  similarities  in 
the  curves  for  the  two  bottom  theories.  One  illustration  is  the  overall  shapes 
of  the  two  solid  curves  for  R/H  =  5.  Also,  we  see  that  except  for  low-density 
bottoms,  Mr  tends  to  increase  as  p2  increases.  This  same  trend  holds  for  the 
MacKenzie  curves.  For  R/H  =  5,  the  peak  Rayleigh  value  in  Fig.  8  is  about 
1 4  dB  higher  at  p2  =  2.10  than  at  p2  =  1.76,  indicating  again  the  generally 
larger  fluctuation  in  received  intensity  for  low-density  bottoms. 

In  Fig.  9,  we  have  taken  R/H  =  20  and  have  graphed  the  limiting  cases  when 
L/H  approaches  zero  and  infinity.  Of  course,  the  MacKenzie  curves  are  the  same 
as  two  of  those  appearing  in  Fig.  3.  There  is  a  substantial  difference  between 
bottoms  with  p2  less  than  about  1 .55  and  those  above  this  value.  For  p2  > 

1  .55,  MacKenzie  results  are  always  greater  than  those  with  Rayleigh,  for  both 
L/H  values;  below  this  value,  the  opposite  is  true.  For  example,  if  L/H  =  0 
and  p2  =  1.45,  the  moment  ratio  for  a  Rayleigh  bottom  is  about  4  dB  greater 
than  that  for  the  MacKenzie  model.  However,  at  p2  =  2.05  and  L/H  =  0,  the 
MacKenzie  MR  is  about  2  dB  more  than  the  corresponding  Rayleigh  value.  In  all 
four  curves,  a  small  correlation  length  L  corresponds  to  less  variance  in 
intensity,  as  discussed  in  conjunction  with  Fig.  3.  Some  bottom-model  consis¬ 
tency  appears  in  Fig.  9,  in  the  sense  that  dB-dif f erence  between  the  two  solid 
curves  (L/H=0)  are  approximately  equal  to  dB-dif ferences  between  the  two  dashed 
curves  (L/H=°°)  at  the  peaks  of  the  Rayleigh  graphs.  For  example,  at  p2  =  1.87, 
there  is  nearly  a  2  dB  difference  between  both  the  L/H  =  0  and  the  L/H  =  « 


curves . 
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The  MacKenzie  results  in  both  Figs-  8  and  9  show  vestiges  of  the  discon¬ 
tinuities  that  appear  in  the  approximate  Rayleigh  curves-  These  appear  as 
local  minima  at  values  of  P2  where  Rayleigh  asymptotes  occur.  Indeed,  if  a/ 8 
were  allowed  to  approach  zero,  the  MacKenzie  minima  would  progressively  deepen. 
Thus,  our  intensity  results  reflect  the  proper  relationship  between  the  two 
bottom  models.  The  MacKenzie  model  is  commonly  regarded  as  the  more  sophisti¬ 
cated  of  the  two  since  it  incorporates  an  additional  physical  effect  (attenua¬ 
tion)  .  Although  certain  relationships  between  the  intensity  results  for  these 
two  bottom  models  have  been  described  and  illustrated,  it  is  also  true  that  the 
levels  of  corresponding  curves  show  some  definite  differences.  For  example, 
the  solid  curves  in  Fig.  9  (corresponding  to  L/H  =  0)  demonstrate  that,  when 
comparing  peak  values  of  the  Rayleigh  curves  with  their  MacKenzie  equivalents, 
there  may  be  up  to  2  dB  differences  between  the  two.  Thus,  the  physical 
features  which  are  properly  incorporated  into  the  bottom  model  can  be  important 
in  determining  the  values  of  the  intensity  moments. 

VI .  SUMMARY 

( 

In  this  paper  we  study  the  ramifications  of  horizontal  random  fluctuations 
in  bottom  structure  on  the  intensity  of  a  received  acoustic  signal.  A  shallow- 
water  sound  channel,  having  constant  sound  speed  and  density,  is  assumed  and 
ray  theory  is  used.  The  water-bottom  interface  is  taken  to  be  a  plane  horizon¬ 
tal  surface,  and  the  ocean  bottom  has  random  density  and  sound-speed  fluctua¬ 
tions  in  the  horizontal  direction  with  a  degree  of  relationship  which  varies 
with  distance.  For  this  stochastic  bottom,  we  derive  expressions  for  ray  path 
length,  spreading  loss,  travel  time,  and  bottom  loss  and  phase  shift.  Expres¬ 
sions  for  the  latter  two  quantities  are  developed  exploiting  perturbation 
expansions,  assuming  small  fluctuations  in  bottom  density  and  sound  speed. 
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Using  these  results,  we  find  the  mean  acoustic  intensity  and  variance  of 
the  intensity  for  a  cw  signal  propagating  in  the  shallow  ocean  channel. 

Stochastic  averaging  over  bottom  randomness  and  incoherent  averaging  over 
per-ray  phases  are  performed  to  obtain  the  intensity  moments.  The  results  are 
sufficiently  general  that  a  variety  of  bottom-acoustic  models  can  be  employed 
with  them.  Also,  the  intensity  moments  are  modeled  to  contain  an  arbitrary 
horizontal  bottom  correlation  coefficient. 

To  illustrate  our  results,  we  chose  bottom-reflection  models  of  MacKenzie 
and  Rayleigh.  In  order  to  use  MacKenzie  theory,  a  procedure  was  devised  to 
specify  values  for  the  attenuation  coefficient.  Also,  our  numerical  results 
assume  a  bottom  coefficient  of  correlation  of  Gaussian  form. 

A  number  of  conclusions  are  drawn  from  our  calculations.  For  example,  we 
find  that  the  standard  deviation  of  received  intensity  is  less  for  sound  trans¬ 
mission  over  fast  (high  density)  bottoms  than  over  slow  (low  density)  bottoms. 
Also,  when  dealing  with  fast  bottoms,  i.e.  mean  densities  over  1.76  g/cm^,  the 
ratio  of  mean  intensity  to  standard  deviation  in  dB  (moment  ratio)  increases  as 
source-receiver  range  R  decreases.  The  same  conclusion  cannot  be  drawn  for 
slow  bottoms,  where  no  distinctive  pattern  exists.  In  general,  we  find  that 
the  moment  ratio  increases  as  mean  bottom  density  increases. 

It  is  shown  that,  due  to  a  cancelling  of  bottom  density  fluctuations,  the 
standard  deviation  of  intensity  increases  as  the  correlation  length  L  increases. 
This  effect  is  a  consequence  of  the  features  of  our  model.  We  find  also  that 
varying  L/H  (H  is  water  depth)  does  not  have  as  great  an  effect  on  the  moment 
ratio  as  does  varying  the  channel  aspect  ratio  R/H. 

We  examine  two  quantities,  normalized  dB-measures  of  intensity  and  inten¬ 
sity  spread,  in  order  to  determine  when  received  intensity  is  within  one 
standard  deviation  of  its  mean  value.  We  find  that  a  greater  spread  of  intensity 
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values  can  be  expected  over  a  low-density  bottom  than  over  a  high-density 
bottom-  Also,  the  spread  in  expected  intensity  values  about  the  mean  increases 
with  R/H. 

We  observe  that,  for  a  given  water  depth  H,  shortening  the  correlation 
length  L  narrows  the  interval  in  which  the  intensity  is  within  one  standard 
deviation  of  its  mean-  It  is  also  shown  that  increasing  the  standard  deviation 
of  bottom  density  fluctuations  causes  a  widening  of  the  spread  of  intensity 
values . 

Statistics  of  intensity  using  the  Rayleigh  model  are  compared  with  those 
found  using  MacKenzie  theory.  The  two  bottom-reflectivity  models  give  results 
that  are  sometimes  similar.  However,  one  of  the  sharpest  differences  is  the 
presence  of  deep  fades  in  moment  ratio,  for  some  mean  bottom  densities,  when 
calculations  are  performed  for  the  Rayleigh  bottom-  These  occur  when  the 
inclination  angle  of  an  acoustic-ray  arrival  in  the  total  field  approaches  the 
critical  angle,  in  which  case  the  variance  of  intensity  becomes  very  large. 

Such  a  behavior,  with  comparable  magnitude,  is  not  observed  when  a  MacKenzie 
bottom  model  is  used. 
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LEGENDS 

(a)  A  horizontal  isospeed  channel  and  the  geometry  of  a  ray. 

(b)  Four  possible  ray  configurations  for  n  =  1  bottom  reflections. 
Moment  ratio  MR  versus  mean  bottom  density  p2  for  four  R/H  ratios: 
L/H  =  1  ,  °z  ~  0*1,  c1  =  1523  m/s,  p-j  =  1  .025  g/cm^.  Surfaced  source 
and  receiver. 

Moment  ratio  MR  versus  mean  bottom  density  p2  for  four  L/H  ratios: 

R/H  =  20.  Other  parameters  as  in  Fig.  2. 

Intensity  ratio  IR  versus  mean  bottom  density  p2:  L/H  —  1,  R/H  =  5. 
Other  parameters  as  in  Fig.  2. 

Intensity  ratio  IR  versus  mean  bottom  density  p2:  =  1,  R/H  «  20. 

Other  parameters  as  in  Fig.  2. 

Intensity  ratio  IR  versus  p2  for  two  L/H  ratios:  R/H  =  20.  Other 
parameters  as  in  Fig.  2. 

Intensity  ratio  IR  versus  p2  for  four  a£  values:  L/H  =  1,  R/H  =  20. 
Other  parameters  as  in  Fig.  2. 

Moment  ratio  MR  versus  mean  bottom  density  p2  using  MacKenzie  theory 
(heavy  curves)  and  Rayleigh  theory  (light  curves)  for  R/H  =  5  (solid) 
and  20  (dashed):  L/H  =  1.  Other  parameters  as  in  Fig.  2. 

Moment  ratio  MR  versus  p2  using  MacKenzie  theory  (heavy  curves)  and 
Rayleigh  theory  (light  curves)  for  L/H  =  0  (solid)  and  00  (dashed): 

R/H  =  20.  Other  parameters  as  in  Fig.  2. 
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